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Single-particle excitation spectra of the two-dimensional Hubbard model on the square lat- 
tice near half filling and at zero temperature are investigated on the basis of the self-consistent 
projection operator method. The method guarantees a high accuracy of the spectra with high 
energy and high momentum resolutions. It takes into account long-range intersite correlations 
as well as the strong on-site correlations. Effects of nonlocal excitations are clarified by com- 
paring the results with those of the single-site approximation. The calculated spectra verify 
the quantum Monte-Carlo results for finite temperatures. The spectra at the Fermi level yield 
a hole-like Fermi surface in the underdoped region and an electron-like Fermi surface in the 
overdoped region. From a numerical analysis of the momentum dependent effective mass and 
self-energy, it is concluded that a marginal Fermi-liquid like state persists even at finite doping 
concentrations in the strongly correlated region because a van Hove singularity is pinned to 
the Fermi surface. It is also found that a kink structure appears in the quasiparticle energy 
band in the same region. The kink is shown to be caused by a mixing between the quasiparticle 
band and an excitation band with strong short-range antiferromagnetic correlations. These re- 
sults suggest an explanation for some of the unusual properties of the normal state in high-Tc 
cuprates. 
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1. Introduction 

The two-dimensional (2D) Hubbard model has been 
widely investigated in the past two decades because it 
has been the simplest model for describing the high- 
temperature superconductivity in layered Cu-based per- 
ovskites.^"^ The Cu 3d^2_j^2 orbitals on a square lattice 
hybridize with the O 2pa; and 2pj^ orbitals on the Cu02 
plane and form an antibonding band near the Fermi level. 
There is a strong on-site Coulomb repulsion of electrons 
in hybridized Cu 3d orbitals with different spins, which 
leads to the 2D Hubbard model. ''"^ 

The 2D Hubbard model at half-filling with nearest- 
neighbor electron hopping is an antiferromagnetic insu- 
lator with the Neel temperature Tn = according to 
the Mermin- Wagner theorem.^ When holes are doped, 
charge fluctuations counteract antiferromagnetism and 
a metallic state should be realized. Rigorous results are 
however scarce especially for intermediate Coulomb in- 
teraction strengths, and therefore the model is of current 
interest in condensed matter physics. 

Numerical techniques applied to the doped regime 
such as the Lanczos method^ and Quantum Monte-Carlo 
(QMC)^"^" have verified the formation of a gap at half- 
filling and clarified global features of excitation spec- 
tra from the insulator to a metal. The antiferromag- 
netic correlations have been shown to reduce rapidly with 
increasing hole concentration. -'^^ More recent calcula- 
tions based on the dynamical cluster approximation^^ 
verified a hole-like Fermi surface in the underdoped re- 
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gion and an electron-like Fermi surface in the overdoped 
region. ^'^ The calculations also suggested the possibilities 
of a pseudogap state as well as the occurrence of super- 
conductivity.^^ 

The methods mentioned above are based on a cluster 
approach. The approach treats a finite cluster explicitly 
by assuming a boundary condition, or implicitly by em- 
bedding it in an effective medium. Since it is numerically 
not possible to treat a large cluster, finite size effects are 
inevitable in these methods. In fact, the methods do not 
yield high momentum and energy resolutions for excita- 
tion spectra, and do not describe the long-range intersite 
correlations due to finite cluster size.^^ It is difficult to 
achieve a momentum resolution better than 0.005 
which was recently obtained in angle resolved photoemis- 
sion spectroscopy (ARPES). This difficulty becomes 
serious for the investigations of the kink structure of 
the quasiparticle band found in the cuprates. Therefore 
one needs different approaches to study a more detailed 
structure of spectra in the 2D Hubbard model. 

In this paper we investigate the single-particle excita- 
tion spectra and related properties of the 2D Hubbard 
model on the square lattice by making use of the self- 
consistent projection operator method (SCPM).^^ The 
SCPM is an extension of the projection operator method 
which uses the coherent potential approximation (PM- 
CPA),^* when determining the momentum dependent 
self-energy. The SCPM treats directly the retarded Green 
function, describes the local on-site correlations by using 
an effective medium, and takes into account the intersite 
correlations up to infinity by using an incremental cluster 
expansion. ^^'^"^ Thus the method yields high momentum 
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and energy resolutions, and allows us to determine the 
detailed structure of nonlocal excitations in both momen- 
tum and energy spaces. Using the SCPM, we will reveal 
the nonlocal effects on the excitation spectra of the 2D 
Hubbard model and its momentum dependent proper- 
ties at zero temperature, which could not be clarified by 
means of the other methods. 

The high-Tc cuprates have unusual properties in the 
normal state; the resistivity shows a T-linear tempera- 
ture dependence and a temperature independent term 
appears in the nuclear relaxation rate.^^ These features 
were first explained by a phenomenological marginal 
Fermi liquid (MFL) scenario.^^ The MFL is defined by 
electrons whose imaginary part of the self-energy near 
the Fermi surface, ImS(a;,T) is proportional to |a;| in- 
stead of up' where w (> T) is an excitation energy and T 
is the temperature. Accordingly, Re T) oc t<;ln|t<;|. In 
spite of successes of the theory, the MFL has remained 
puzzling because no microscopic theory could be pro- 
vided for it. We show in the present paper that a MFL 
like behavior appears in the underdoped region, which 
may be considered as a form of justification of the phe- 
nomenological MFL theory. 

Another interesting point of our results is that a kink 
structure in the quasiparticle band appears in the under- 
doped region. Recent ARPES experiments have revealed 
a kink in the quasiparticle band of higli-Tc cuprates in 
both the normal and superconducting state. ^^"^^ The 
typical kink energy at which the quasiparticle velocity 
changes is 70 meV. Previous theories for the kink struc- 
ture assumed either a coupling of electrons to a magnetic 
resonance mode found in inelastic neutron scattering ex- 
periments^^' or a coupling to a longitudinal optical 
phonon mode.^^'^^ These theories do not take into ac- 
count other electron correlations which are an essential 
ingredient of the high-Tc cuprates. We demonstrate on 
the basis of the SCPM that the kink structure can be 
explained by electron correlations without introducing 
lattice degrees of freedom. This provides a new possible 
mechanism for the observed kink structure in cuprates. 
Preliminary results on these topics have been published 
recently. ^^'^^ It should be noted that the validity of nu- 
merical calculations on these problems strongly depends 
on the momentum and energy resolutions. We therefore 
have recalculated the spectra increasing the number of 
mesh from 80 x 80 to 160 x 160 in the first Brillouin zone 
and the energy mesh by a factor of two. On the basis 
of the new results of calculations we will discuss in the 
present paper the MFL and the kink structure in more 
details. 

The paper is organized as follows. In the next section, 
we review briefly the SCPM to the Hubbard model. In 
Sec. HI we present our results of calculations for the ex- 
citation spectra. We limit ourselves here to excitation 
spectra in the normal state at zero temperature. In Sec. 
HI A we discuss the characteristics of the excitation spec- 
tra at half-filling under the assumption of a paramagnetic 
effective medium. In Sec. HI B, we treat the doped case. 
We present the excitation spectra, integrated density of 
states (DOS), Fermi surface, momentum distribution, 
and the momentum dependent effective mass. We clar- 



ify the nonlocal effects of electron correlations on these 
quantities by comparing them with those of the single- 
site approximation (SSA) and other numerical results. 
The MFL-like behavior which we find in the underdoped 
region is discussed separately in Sec. HI C. The MFL 
features are shown to be caused by a pinning of the van 
Hove singularity to the Fermi surface under doping. In 
Sec. HI D, we discuss the kink structure of the quasipar- 
ticle band, which is obtained in the underdoped region. 
The results of our calculations are summarized in Sec 
IV, where also the unusual behavior of cuprates in the 
normal state is discussed. 

2. Self-consistent projection operator method 

We consider a Hubbard model on a square lattice with 
an atomic level eo; a nearest-neighbor transfer integral 
tij (= t), and an on-site Coulomb interaction parameter 
V. 



H 



(1) 



Here the chemical potential fi has been added to the 
Hamiltonian. Furthermore {dia) is the creation (an- 
nihilation) operator for an electron with spin a on site i 

and Uia = alaO-ia- 

The single-particle excitation spectrum is obtained 
from the poles of the retarded Green function given by 
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Here z = Lu + iS with S being an infinitesimal positive 
number, Cka is the Hartree-Fock energy given by eka- = 
eo + ek — H + U {rii-cr) where {ni^) and Ck denote the 
average number of electrons with spin a on site i and 
the dispersive part of the energy band for noninteracting 
electrons, respectively. 

The self-energy Akcr{z) is given by a Fourier transform 
of the memory function Mij^{z) as 



^ka{z) = U'^Y1 ^jOa{z) exp{ik ■ Rj) . 
j 



(3) 



In the projection operator method,^*^ the memory 
function is given by 



aW=(4. I {z-L)-'A]^) . 



(4) 



j)). The inner product of op- 



Here = 

erators A and B is defined by (A|_B) = 5] + ) where 

the expectation value is taken with respect to the ground 
state of the system. The dynamics of the electronic sys- 
tem is described by a superoperator L. It acts on other 
operators A according to LA = [H, A\ _ . Furthermore 
L = QLQ where the projector Q is given by Q = 1 — P 
and P = Y.ia l«L)('^Ll- "^^^^ Q eliminates \a\^) from 
further considerations. 

In the SCPM^^ we introduce an effective Liouville op- 
erator 



L{z)A=[H{z),A]. 



H{z) = Ho 



(5) 
(6) 
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Here Hq is the Hartree-Fock Hamiltonian, and 'Scr{z) de- 
fines an energy-dependent effective medium and is called 
the coherent potential. The Liouvillean L is divided into 
a coherent part L{z) and an interaction part Li{z). By 
making use of the multiple scattering theory, we can de- 
rive an incremental cluster expansion for the memory 
functions: 



+ 1^ ^ AM<>")(.) + ... , (7) 



+ 



1 



with 



AMf\z)=Mt^iz) 

-AM^{z)-I^mI^{z)-M^{z) 



AM^{z) = M^'{z)-M]^{z) 



im). 
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(9) 

(10) 
(11) 



AM^-^(z) 



M^-^(z) 



-AM^fiz) - AMgrH^) - M^;i>{z) . (12) 
In the above expressions, M^'^l.{z) are cluster memory 



functions defined by 



M, 
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Here l''''\z) = QL^°\z)Q, and Sriia- = Uia- - {riia). 
L'^'^\z) is the Liouvillean for the cluster c which is cm- 
bedded in an effective medium Yi„{z)^ and denotes 
the number of atoms in the cluster. It should be noted 
that we take into account all the 'clusters' obtained by 
choosing any sites from the 2D square lattice points. 
This is quite different from the usual cluster approaches 
in which only one compact cluster of sites with a given 
form is considered. It should also be noted that when 
the off-diagonal memory function (eq.(8)) is neglected 
and only the single-site term (SSA) is taken into account 
in eq. (7), the SCPM reduces to the PM-CPA, which 
is equivalent to the many-body CPA,^^ the dynamical 
CPA,^^ and the dynamical mean-field theory. ■^^''^^ 

We use here the memory functions obtained by the 
renormalized perturbation scheme: 



J ija 



(15) 



The diagonal matrix h'^f^ is defined by^^ 

^(i) ?7(1 - 2(n,_,)) 



(16) 
(17) 



The screened memory function Mq in eq. (15) is a A^c x 
Nc matrix given by 



M^liz) 



d^d^d^'P&i^P^.i^)f^.(^')x{e, e\ e") ^ ^^^^ 



e' + e" 



(19) 



(ni-o-)o(l - {ni-a)o) 
X{e,e',e") = f{-e)f{-e')fie") + f{e)fie')fi-e") . (20) 

^ (c) 

The memory function Mg^j^ consists of that of the 
second-order perturbation^^' theory and of a prefac- 
tor Aij^ ensuring the correct second moment in the mo- 
ment expansion. Furthermore {ni^)^ is the average elec- 
tron number defined by {nia-)o = / d'^f{'^)Pul-{'^) where 
/(e) is the Fermi distribution function. The density of 

states Pijlie) is defined by 



{F{z)) 
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(22) 



Co- - Scr(z) - e 

where Ca- = cq — IJ- + U (rii-o-) , and pij (e) is the density of 
states for a noninteracting system, 

Pij{e) = -^Yd{e- ek)exp[-ik- {Ri- Rj)] . (23) 



k 



The coherent potential introduced in eq. (6) is deter- 
mined self-consistently from the CPA equation, i.e., 



t^z) = Aua{z) = N-'Y.^kaiz) . 

fc 



(24) 



Equations (3), (7), (8), (15), and (24) form a self- 
consistent set of equations from which we can obtain 
the effective medium ^^(z). The nonlocal excitations 
are then obtained from eq. (2). We have solved the self- 
consistent equations at zero temperature within the two- 
site approximation in which we take into account the first 
two terms of the right hand side of eq. (7) and the first 
term in eq. (8). 

3. Numerical results 

3. 1 Nonlocal excitations of the half-filled band 

The 2D Hubbard model is considered to be an antifcr- 
romagnetic insulator with the Neel temperature 7n — 
K at half filling.^ In the effective medium approach there 
are two solutions found, a nonmagnetic solution and an 
antiferromagnetic (AF) one. The latter describes the sys- 
tem at half filling. But in the SCPM an AF effective 
medium overestimates antiferromagnetic correlations at 



4 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



0.3 




Energy 



Fig. 1. Singlc-particlc excitation spectrum along the high sym- 
metry lines for half filling (n = 1) and U = S. The Fermi level is 
indicated by a bold line. Here and in the foUowings, the energy 
is measured in unit of |t|. Note that F = (0, 0), X = (tt, 0), M = 
{it, tt) are in unit of the lattice constant. 




Fig. 2. Single-particle excitation spectra in the single-site approx- 
imation (SSA). Parameters are the same as in Fig. 1. 

finite doping concentrations because the medium elimi- 
nates spin fluctuations. Therefore the SCPM does not de- 
scribe properly the AF at T = 0. Here we study the prop- 
erties of excitations in the presence of the nonmagnetic 
medium, because the nonmagnetic solution at half-filling 
is smoothly connected to that at finite doping concentra- 
tion. In the following we present results for the nonlocal 
excitations in the nonmagnetic medium at half-filling. 

For the numerical calculations, we adopted a 160 x 160 
(or 80 X 80) mesh in the first Brillouin zone, and obtained 
the momentum dependent self-energy (3) by taking into 
account the off-diagonal memory functions up to the 50th 
(!) nearest neighbors. The latter are calculated from eq. 
(15). The screened memory functions MQ^j^(2) in eq. (15) 
are calculated by means of the Laplace transformation 
method. In order to avoid singularities in the Green 
functions, we used a complex energy z = uj + iS with 



Fig. 3. Total densities of states (DOS) at half filling as function 
of U. Noninteracting DOS is shown by dotted line. 

S = 0.025 (or 6 = 0.05) in unit of \t\. 

Figure 1 shows the excitation spectrum along the high 
symmetry line at U = 8. The spectrum is characterized 
by incoherent Mott-Hubbard excitations, i.e., a lower 
Hubbard band (LHB) around the F point and an upper 
Hubbard band (UHB) around the M point. In addition 
there is the quasiparticle band near the Fermi energy 
{to = 0). The nonlocal excitation spectrum is compared 
with that in the single-site approximations (see Fig. 2). 
We find that nonlocal correlations reduce the amplitude 
of the quasiparticle peak near the Fermi level, and in- 
crease its band width by a factor of two as compared with 
the SSA for an intermediate strength of the repulsion U. 
Alternatively, the electron correlations in the SSA over- 
estimate the band narrowing in the quasiparticle state. 
When intersite correlations are additionally included, the 
electrons recover parts of the original degrees of freedom 
for their motion. The result is a band broadening but the 
Coulomb repulsion energy remains optimized. 

Nonlocal excitations also enhance the amplitude of the 
LHB around the F point as well as that of the UHB 
around the M point. Similar global features of these non- 
local excitation spectra are also found in the 3D Hubbard 
model on a simple-cubic lattice.^'' Another characteris- 
tic of the spectra is pronounced flat-band excitations at 
LV « ±2.0 due to antiferromagnetic correlations. These 
flat bands are precursors of an antiferromagnetic insu- 
lator.^^ Thoy should develop further when higher-order 
cluster correlations are taken into account by going be- 
yond the two-site approximation. 

The total density of states (DOS) as a ftmction of in- 
teraction U is presented in Fig. 3. A characteristic fea- 
ture of the DOS is a logarithmic divergence at a; = 0, 
which is due to a van Hove singularity at the X point. 
For a small Coulomb interaction, we can compare our 
results with those of the perturbation theory. We 
find a good agreement between the presend DOS and 
that obtained by the 4th order perturbation theory for 
U = 2. For U = 4, the shoulder position is |a;| = 5.0 in 
the present calculation, while it is |ci;| = 4.5 in the 4th or- 
der perturbation. These results verify the validity of our 
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Fig. 4. Sclf-cncrgics at the Fermi momentum (7r/2,7r/2) in var- 
ious approximations; Two-site approximation (solid curves), 
second-order perturbation (thick dashed curves), single-site ap- 
proximation (thin dashed curves). Note that the real part is an- 
tisymmetric, and the imaginary part is symmetric in energy. 

results for the regime of Coulomb interaction W^A, in 
which the 4th order perturbation is believed to be valid. 
The quasiparticle band in the low energy region becomes 
narrow when the interaction U is increased. The Mott- 
Hubbard incoherent bands on the other hand develop in 
the high energy regions. Furthermore, the nonlocal anti- 
ferromagnetic correlations create subbands at » 1.8 
when ;7>6. When C/>10, we find within the two-site 
approximation a negative spectral density in a small re- 
gion of momentum space. Therefore we have to include 
higher-order cluster correlations in order to discuss the 
spectra beyond within the present scheme. In the 

following we limit our discussions to the metallic region 
U<10. 

A characteristic fcattirc of the 2D Hubbard model at 
half filling is that of marginal Fermi liquid. The MFL is 
defined by Re Afe^(z) oc wlnjo;! and Im Akp{z) oc \uj\ 
for small and for the Fermi wave number kp- In 
the weakly correlated regime such behavior is found by 
2nd-order perturbation theory.'^*' "^^ For the intermedi- 
ate Coulomb interaction strength, we have calculated the 
self-energy at the Fermi surface. An example is shown in 
Fig. 4. The MFL behavior seems to be for k present even 
for strong Coulomb interaction U, although numerically 
the linear dependence of Im Ahp{z) on is not so ac- 
curate. 

We have calculated the momentum dependent effec- 
tive mass nik = (1 — dIieAk{z) / duj)i^=o using numerical 
differentiation. As presented in Fig. 5, the calculated ef- 
fective mass shows a strong momentum dependence, hav- 
ing a minimum at the M point (= (tTjTt)), and a maxi- 
mum at the X point (= (7r,0)). Along the Fermi surface 
cos kx+cos ky = 0, nik has a minimum at {n/2, ■jt/2) and 
a maximum at (tt, 0). This suggests the appearance of 
the 'Fermi arc'^^ when holes arc doped. 

Because ReAfe(2;) oc a;ln|cL-'| for small \u!\ in the MFL, 
it is irik — 1 oc —ln6u} for small energy steps Suj when a 
numerical differentiation is done. Therefore, the simple 
numerical method for the calculation of fc-dependent ef- 




Fig. 5. Momentum dependent effective mass at half filling, calcu- 
lated by numerical difi^erentiation by using an energy interval of 
Sui = 0.05. 




Fig. 6. Excitation spectrum along the high symmetry lines for 
doping concentration Sf^ = 0.14. 

fective mass breaks down in the case of the MFL. We 
have verified that the effective mass at the Fermi surface 
in Fig. 5 increases logarithmically with decreasing Suj, 
i.e., m/c along the Fermi surface is expected to show a 
logarithmic divergence. We will return to this problem 
in Sec. Ill C. 

3.2 Nonlocal excitations in the doped Hubbard model 

When holes are doped into a half-filled 2D Hubbard 
system, quantum spin and charge fluctuations are en- 
hanced due to a motion of holes, and thus the system is 
expected to lose the long range magnetic order. In the 
very weak Coulomb interaction limit, it is suggested by 
using the Hartree-Fock approximation^^ that an incom- 
mensulate antiferromagnetic order can remain even at 
finite doping concentration, and the insulating state can 
survive under the long range magnetic order. A problem 
of the theory is that a finite value of Tn at half-filling con- 
tradicts with the Mermin- Wagner theorem.^ This sug- 
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Fig. 7. Contour map of the excitations at (5^ = 0.05. Closed cir- 
cles with error bars are the QMC results^*^ at T = 0.33. The 
dashed curves show the Hartree-Fock quasiparticle dispersion. 



Fig. 8. Contour map of the excitations at <5^ = 0.14. 



15 




gests that the quantum fluctuations are dominant in the 
2D system so that the long range magnetic order is sup- 
pressed at finite doping concentration. Although there 
are no rigorous results on this problem we assume here 
in our numerical calculations the nonmagnetic state even 
for an infinitesimal doping concentration. 

A typical excitation spectrum is shown in Fig. 6 for 
a doping concentration of Sh = 0.14. As seen from the 
figure, the upper Hubbard band moves upwards. The 
lower Hubbard band has become weaker and merges with 
the quasiparticle band near the T point. The excitations 
due to short-range AF correlations, which are located at 
« 1.8 in Fig. 1, have disappeared. Consequently, the 
spectral weight of quasiparticle peak is enhanced. 

We have examined the doping dependence of the ex- 
citation spectra. At 5 % doping, the LHB is nearly de- 
stroyed as shown in Fig. 7. The flat bands at |a;| ~ 2 
around the F and M point, which are due to short-range 
antiferromagnetic correlations, have been weakened. The 
quasiparticle band on the other hand is well developed. 
Accordingly, the flat quasiparticle band around the X 
point sinks below the Fermi level. In the optimum doped 
region (see Fig. 8), the quasiparticle band is further pro- 
nounced. But its width is broadened by 26 % as com- 
pared with 5 % doping. Note that due to the increase of 
holes the flat band around the X point is at the Fermi 
level. In the overdoped region (see Fig. 9), the quasipar- 
ticle band gains further weight and its width increases 
by 10 % as compared with 14 % doping. The flat band 
around the X point is now above the Fermi level. 

We have compared our results at zero temperature 
with those of the QMC at flnite temperatures.^*^ As 
shown by closed circles in Figs. 7, 8, and 9, the present re- 
sults for the nonlocal excitations are consistent with the 
QMC results for the underdoped region as well as for 
the overdoped one. Especially, the calculated quasipar- 
ticle bands show quantitative agreement in both cases. 
In the QMC calculations for 20 % doping concentration, 
weak excitations are found at w w —3.0 around the X 
and (7r/2,7r/2) points as shown in Fig. 9. Although not 
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Fig. 9. Contour map of excitation spectra at Sh = 0.20. 



visible in Fig. 9, we have verified that the same exci- 
tations do appear in our calculations. For example, we 
find a small peak of height 0.03 at the X point (tt, 0) 
and energy lu = —2.9, in agreement with the QMC re- 
sults. We have also investigated the spectra in the SSA. 
The quasiparticle weight in the SSA is underestimated 
at the F point and is overestimated at the M point. The 
LHB (UHB) in the SSA is correspondingly excessively 
enhanced (weakened). 

The integrated DOS's are shown in Fig. 10 for differ- 
ent hole dopings 5h = 0.05, 0.14, and 0.20. We find that 
the weight of the LHB is rather small for these doping 
concentrations. The upper Hubbard band simply shifts 
to the higher energy region with increasing doping con- 
centration. The sharp quasiparticle peak near the Fermi 
level is well developed. The peak first sinks and gradually 
rises with increasing doping concentration. This behav- 
ior is directly connected with the change of the flat band 
near the X point with hole doping. 

We have also investigated the excitation spectrum at 
w = 0. Results are shown in Figs. 11 and 12. We flnd a 
hole-like Fermi surface in the underdoped region and an 
electron-like Fermi surface in the overdoped region. This 
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Fig. 10. Total DOS for various doping concentrations 

Fig. 12. Excitation spectrum at the Fermi energy for Sf^ = 0.20. 



p(k,a>=0) 




5h=0.00 
8h=0.05 
8h=0.14 
8h=0.20 
Sh=0.25 




Fig. 11. Excitation spectrum at the Fermi energy for <5(, = 0.05. 

behavior is in agreement with the one obtained by the dy- 
namical cluster approximation at finite temperatures.^'^ 
It should be noted that the present results do not satisfy 
Luttinger's theorem. ^'^ The deviation from the volume 
predicted by that theorem is small in the overdoped re- 
gion, but it is large in the underdoped region. The ratio 
of the excess volume to the predicted one is 5v/vq = 0.3 
at 6h = 0.05, for example. It is c;omparable to the value 
of 0.4 obtained by the QMC calculations.^" 

The momentum distributions along the high sym- 
metry lines are presented in Fig. 13. At the M point Uk 
decreases with hole doping as expected. The momentum 
distribution at the V point, on the other hand, decreases 
first and then gradually increases with hole doping. The 
reduction in the underdoped region originates from the 
reduction of the LHB as seen in Fig. 10. It also yields an 
enhancement of at the X point. It should be noted 
that there is no jump of rife at the X point for half- filling 
because the quasiparticle weight vanishes there as will be 
discussed in the next section. It is also interesting that no 
discontinuity of rife is seen at the X point for 5h = 0.14 at 
which the Fermi surface is not well defined because the 



Fig. 13. Momentum distribution along the high symmetry direc- 
tions for various doping concentrations. 

flat band is here right at the Fermi level. These anomalies 
are not found along the F-M line; we find there a clear 
discontinuity of rife at the Fermi surface irrespective of 
the doping concentration. 

3.3 Marginal Fermi liquid behavior 

We have discussed in Sec. Ill A the strong k depen- 
dence of the effective mass and we have pointed out 
its logarithmic divergence at the Fermi level for the case 
of half filling (see Fig. 5). When holes are doped, the ef- 
fective mass rapidly looses the k dependence. We show in 
Fig. 14 an example for 5 % doping, where the maximum 
of TOfe changes from the X point to the F point. The 
maximum and minimum values of rrife with increasing 
doping concentration are shown in Fig. 15. Both quan- 
tities increase with decreasing doping concentration. At 
(5^ = 0.025, we find a sudden change of m^. When the en- 
ergy in the numerical differentiation is changed in steps 
varying from 5u} = 0.05 to 5uj = 0.005, the maximum 
value of rrife at the X point increases logarithmically for 
6h < as shown in Fig. 15. This implies that up to the 
finite doping concentrations of 2.5 % the electrons are 
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Fig. 16. Self-energy at the Fermi momentum k = (7r,0.137r) for 
doping concentrations Sh = 0.02 and 0.05. 
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Fig. 15. Momentum-dependent effective mass vs. doping con- 
centration (5^. Open circles: maximum value at the X point for 
Sh ^ 0.02 and F point for ^ 0.02, closed circles; minimum 
value at the M point. Numerical derivatives are taken with en- 
ergy steps of 6(jj = 0.05. In the case of maximum m^, the results 
for Su! = 0.005 are shown by +. The momentum-independent 
effective mass in the single-site approximation (SSA) is shown 
by the dashed curve. 



forming a MFL. We can not verify numerically though 
how accurately a MFL is realized at those concentra- 
tions. Figure 16 shows the calculated self-energies at a 
particular kp for 6h = 0.02 (in the MFL region) and 
Sh = 0.05 (in the Fermi liquid region). The Im Akp{z) at 
6h = 0.02 is approximately proportional to |a;| for small 
LO. We therefore conclude that the quasiparticle states in 
the low doping region arc at least very close to a MFL. 
We note that a MFL can not be described by the conven- 
tional cluster approaches because the long-range intersite 
correlations near the Fermi surface are indispensable for 
it to occur and they are not taken into account in the 
usual cluster theories. 

As mentioned before, the change from the Fermi liq- 
uid (FL) to the MFL takes place discontinuously. This 
is clearly seen in the doping dependence of the chemi- 



cal potential shown in Fig. 17. The chemical potential 
increases monotonically with decreasing doping concen- 
tration and jumps at = 0.025 by A/x ~ 0.55. Between 
(5/i = 0.020 and 6h = 0.025, we find two solutions. Below 
6h = 0.020, there is only one solution corresponding to a 
MFL. 

The discontinuity in the chemical potential as well as 
in TOfe at 5^ originates from the collapse of the LHB 
with increasing hole doping. Figure 18 shows the total 
DOS at Sh = 0.02 and 0.05. At low doping concentra- 
tions, the spectral weight of the LHB moves directly to 
the UHB with increasing doping because doubly occu- 
pied states are almost totally suppressed when correla- 
tions are strong. Consequently, the van Hove singularity 
is pinned to the Fermi level as shown in Fig. 18. It leads 
to the MFL behavior at finite doping concentrations up 
to 2.5 %. Further doping, however, causes the LHB to 
disappear and so the peak at the Fermi level, which is 
associated with a flat band around the X point, is moving 
down below the Fermi level. Finally, the spectral weight 
of the LHB moves to the quasiparticle states, and the 
electrons obtain more and more itinerant character with 
reduced correlations. 

Calculated discontinuity lines in the U — Sh plane are 
presented in Fig. 19. The MFL-like state with the well- 
defined LHB is located around Sh = 0.02 and [7 = 8. 
The MFL is separated from the normal FL state with 
collapsed LHB by a hatched region in which two self- 
consistent solutions exist. Note that the hatched region 
is slightly reduced as compared with our preliminary re- 
sult^* because we have increased the numerical acxniracy. 
We also found two additional discontinuity lines in the 
U — Sh plane. One is around Sh = 0.005 and U = 7.5, 
while the other is around Sh — 0.10 and U = 6.5. The 
discontinuities are small (A/x ~ 0.05) at both lines and 
the latter changes to the kink-like anomalies which are 
indicated in Fig. 19 by the dotted lines. We have exam- 
ined these discontinuities of A/x by changing the mesh 
of ui and k in the numerical calculations. Our conclu- 
sion is that the discontinuity with the hatched region and 
that around {Sh, U) = (0.005, 7.5) exists. But, for the line 
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Fig. 18. The DOS of the marginal Fermi liquid like state (solid 
curve) and the normal Fermi liquid state (dot-dashed curve). 
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Fig. 19. Phase diagram showing the marginal Fermi liquid regime 
(MFL) and the normal Fermi liquid state (Normal FL) regime. 
The lines at which the chemical potential changes discontin- 
uously arc shown by solid curves. The region with two self- 
consistent solutions is shown by the hatched area. Dashed lines 
along the hatched area arc extrapolations. Along the dotted lines 
the chemical potential has a kink. 




Fig. 21. The excitation spectrum along the nodal direction for 
Sh = 0.05. 



around {6h,U) = (0.10,6.5), the value A/x and position 
seem still to depend on the number of mesh points even 
when both the momentum steps and the energy steps are 
decreased down to n/8Q and 0.025, respectively. There- 
fore we do not exclude the possibility that the line disap- 
pears when the numerical accuracy is increased further. 



3.4 Kink structure in the underdoped region 

We have examined the nonlocal excitations in the MFL 
region. The excitations in this region are similar to the 
ones of the half-filled case (sec Fig. 1). The LHB excita- 
tions are located around the F point at w « —6.0. We also 
find flat-band excitations at a; = ±2.0, which are due to 
short-range antifcrromagnetic correlations. The detailed 
structure of the quasiparticle states is however unusual. 
Figure 20 shows a spectrum along the nodal direction 
for Sh = 0.01. Notice that a kink structure appears in 
the quasiparticle states at (|A;|, w) = (0.527r, —0.85). The 
kink is caused by a mixing of quasiparticle states with 
the excitations due to magnetic short-range order. The 
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Fig. 22. The excitation spectrum along the antinodal direction 
for 5h = 0.01. 



latter are enhanced by the nesting of the Fermi surface in 
the MFL region. When 6h, > 0.03, the LHB collapses as 
mentioned before and the flat band around the X point 
sinks. Accordingly the excitations due to antifcrromag- 
netic correlations are weakened. The quasiparticle state 
is then extended to the higher energy region and the kink 
structure disappears as shown in Fig. 21. 

We have also examined the quasiparticle excitations 
along the antinodal direction. As seen in Fig. 22 the mix- 
ing of the quasiparticle band with the excitations due 
to the short-range antiferromagnetic correlations is no 
longer so strong, and we do not find a clear kink along 
this direction. 

The ratios of the Fermi velocity above the kink energy 
Wkink to that below the one are vp'/vp = 2.7 {5h = 
0.0), 1.8 {5h = 0.01), and 1.5 {Sh = 0.02), respec- 
tively. The concentration dependence of the Fermi ve- 
locity above Wkink is shown in Fig. 23. Along the F-M 
direction, the Fermi velocity is reduced by a factor of two 
as compared with the one of the noninteracting system. 
Its concentration dependence is rather weak. The Fermi 
velocity along the F-X-M line is strongly influenced by 
the energy position of the flat band at the X point. It 
first increases with the shift of the Fermi momentum fc^ 
towards the M point, then decreases gradually with the 
return of kp to the X point. When 5^ = 0.14, kp crosses 
the X point, and the Fermi velocity starts to increase 
again with hole doping. 

4. Summary and discussions 

Wc have investigated numerically the single-particle 
excitation spectra of the 2D Hubbard model at zero 
temperature and various doping concentrations, thereby 
using the self-consistent projection operator method 
(SCPM). The SCPM describes the nonlocal excitations 
due to long-range intersite correlations, and allows us to 
compute the spectra with high energy and momentum 
resolution. The method therefore is able to clarify de- 
tails of nonlocal excitations that had not been resolved 
by previous approaches. 

We have calculated the momentum-dependent excita- 
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noninteracting system are shown by the dot-dashed curve and 
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tion spectra from the underdoped to the overdoped re- 
gion within the two-site approximation. Pairs of sites, 
up to the 50th nearest neighbor were considered. The 
spectra are characterized by a LHB around the F point, 
and an UHB around the M point. They show strongly 
renormalizcd quasiparticle states near the Fermi level, as 
well as excitations caused by magnetic short-range cor- 
relations. We verified that the present results are consis- 
tent with those obtained by QMC calculations at finite 
temperatures. ^° Especially the quasiparticle excitations 
show a quantitative agreement with the QMC results. We 
have shown that compared with the single-site approx- 
imation intersite correlations suppress the quasiparticle 
weight in the vicinity of the Fermi level by a factor of two. 
They relax also the band narrowing of quasiparticles by 
a factor of two for an intermediate strength of the repul- 
sive interaction. Intersite correlations also produce exci- 
tations due to short-range antiferromagnetic order in the 
underdoped region. Around the X point the correlations 
smoothen a sharp drop of the momentum distribution rik 
at the Fermi surface. 

The antibonding band of Cu 3da,2_j^2 orbitals hybridiz- 
ing with the planar O 2pa; and 2py orbitals is approxi- 
mately described by a 2D Hubbard model.^'^'^^ There- 
fore the present results suggest an explanation of some 
of the unusual properties in the high Tc cuprates. We ob- 
tained a hole-like Fermi surface for the underdoped re- 
gion and an electron-like Fermi surface for the overdoped 
region, being in agreement with the experimental data 
on LSCO.'*"'' We found numerically a MFL-likc behavior 
at finite doping concentrations by analyzing the momen- 
tum dependence of the effective mass. The results justify 
a phcnomcnological MFL description^^ which explains 
many aspects of the normal state in cuprates. We found 
that the MFL is caused by a pinning of the van Hove 
singularity to the Fermi level. It is due to a transfer of 
spectral weight from the LHB to the UHB for small dop- 
ing concentrations when U/\t\ is sufficiently large. The 
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MFL region calculated here is limited to a small area 
in the U/\t\ vs. 5h plane. It remains work for the future 
to extend the MFL region by going beyond the two-site 
approximation and by improving the model. 

The present theory yields a large deviation from the 
Luttinger theorem^^ in the underdoped region in agree- 
ment with previous investigations."'^''' -'^^ In our calcula- 
tions, the violation is caused by a sudden reduction of the 
LHB with increasing doping concentration or Coulomb 
interaction strength U . The collapse of the LHB strongly 
pulls down the flat band around the X point and in- 
creases considerably the volume below the Fermi sur- 
face. Accordingly, the spectral weight of the LHB moves 
to a higher energy region. It has not yet been clari- 
fied however, under which general conditions the Lut- 
tinger theorem can be violated. Experimentally, LSCO 
has been reported to satisfy Luttinger's theorem,^^ while 
Ca2-aNa.cCu02Cl2 was shown to violate it strongly .^^'^^ 
In the latter case, it has been observed that the LHB ap- 
proaches the Fermi level with hole doping. This might 
correspond to the collapse of the LHB and the transfer 
of spectral weight to the higher energy region found in 
our calculations. 

We have also found in the underdoped region a kink in 
the excitation spectrum along the nodal direction. The 
kink is caused by a mixing between the quasiparticle 
band and the excitation band caused by short-range an- 
tiferromagnetic correlations. This suggests that the kink 
found in the cuprate is of electronic origin; in principle 
there is no need to include any other (e.g. phononic) de- 
grees of freedom coupled to electrons in order to explain 
the kink structure. In fact, we obtain kink energies of 
'^kink = 130-170 meV when we adopt a reasonable trans- 
fer integral \t\ = 0.2 eV.^^ They are comparable to the 
experimental values of 60-70 meV.^^ Calculated ratios 
of the Fermi velocity vf'/vf at Wkink are between 1.5 
and 2.7 in the underdoped region. These values should 
be compared with the experimental ones,^'' i.e., 1.3-3.0. 
We suggest that a strong doping dependence of the kink 
in La2-2;Sr2;Cu04 (0 < .t < 0.07)''^ can be explained 
by a strong doping dependence of the antiferromagnetic 
short-range order as found in the present calculations. 

The 2D Hubbard model considered here is the simplest 
model for cuprates. It is leading to a Neel temperature 
Tn = at half filling. In order to describe the layered 
Cu-based perovskites with high we have to adopt 
more realistic models, i.e., we have to include 3D fea- 
tures. Furthermore the present calculations neglect the 
self-consistency of the off-diagonal matrix elements of the 
effective medium. Because of the missing self-consistency, 
the results presented here may underestimate the antifer- 
romagnetic correlations. The self-consistency of the off- 
diagonal matrix elements of the effective medium and the 
higher-order cluster correlations should extend the MFL 
and the kink region on the phase diagram. Correspond- 
ing calculations are left for future investigations. 

Note added in proof — Recently a high energy anomaly 
called "waterfall" has been found at about 0.3-0.5 eV 
in the ARPES of high-Tg cuprates [ W. Meevasana et. 



al., cond-mat/0612541vl; J. Graf et. at, Phys. Rev. 
Lett., 98, (2007) 067004; D.S. Inosov et. al, cond- 
mat/0703223vl ]. According to the data on LSCO, the 
low-energy kink anomaly at 60-70 meV strongly de- 
pends on the doping concentration, while the high en- 
ergy anomaly does not. We suggest that the former cor- 
responds to the kink anomaly we found here because of 
the strong concentration dependence and that the latter 
might be a phenomenon related to a mixing between the 
quasiparticle band and the incoherent LHB. 
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